Evidence of 13 > from global neutrino data analysis 
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The neutrino mixing angle #13 is at the focus of current neutrino research. From a global analysis of 
the available oscillation data in a 3u framework, we previously reported [Phys. Rev. Lett. 101, 141801 
(2008)] hints in favor of #13 > at the 90% C.L. Such hints are consistent with the recent indications 
of Vfj, — > v e appearance in the T2K and MINOS long-baseline accelerator experiments. Our global 
analysis of all the available data currently provides > 3u evidence for nonzero 613, with la ranges 
sin 2 #13 = 0.021 ± 0.007 or 0.025 ± 0.007, depending on reactor neutrino flux systematics. Updated 
ranges are also reported for the other 3v oscillation parameters (5m 2 , sin 2 612) and (Am 2 , sin 2 623). 

<n : 

OJO' I. INTRODUCTION 

3 : 

Neutrino oscillation experiments have established that neutrino flavor and mass states do mix In the simplest 
' framework, the three flavor states (f e , t/^, i/ r ) are quantum superpositions of only three light mass states (v\, 1/3) 
£N| , via a unitary mixing matrix U, parametrized in terms of three rotation angles (#12, #13, $23) and one possible CP- 
violating phase 6 in standard notation [l[. The amplitudes and frequencies of flavor oscillation phenomena are 
I governed, respectively, by the 8ij angles and by two squared mass differences, namely, 

Sm 2 = m\ - m\ > (1) 

and, in our convention 

where Am 2 > (< 0) corresponds to normal (inverted) mass spectrum hierarchy. Typically, a single experiment is 
qq ' mainly sensitive to only one of the above mass gaps and to one mixing parameter, although sublcading effects driven 
£S) , by the remaining parameters may become relevant in precision oscillation searches. 

So far, solar and long-baseline reactor neutrino experiments have measured the mass- mixing parameters (5m 2 , 612) 
^sO in the v e — > v e channel, while atmospheric and long-baseline accelerator (LBL) experiments have measured (Am 2 , #23) 
\Q . in the — > channel. Conversely, short-baseline reactor experiments, mainly sensitive to (Am 2 , #i3),have set 
upper but not lower — bounds on the mixing angle #13; see [l[ for an overview. However, we observed in [H, H[ that 
the two data sets mainly sensitive to <5m 2 and to Am 2 provided two separate hints in favor of #13 > which, in 
combination, disfavored the null hypothesis #13 = at 90% C.L. 

The statistical significance of the hints, as well as their possible origin in subleading oscillation effects, have been 
examined in detail and also debated in a number of analyses @-[l3], often triggered by new input data and, more 
recently, also by a new, critical evaluation of older inputs for the reactor neutrino fluxes [HI (see [l9l - l22j ). Within 
the standard 2>v framework (with no extra sterile neutrinos i> s ), the overall statistical significance of #13 > has, so 
far, not exceeded the level of ~ 2er, with a corresponding estimated range sin 2 613 ~ 0.02 ± 0.01, see e.g. (TlLl22T|. 

Very recently (June 2011), new relevant results have been announced by two long-baseline accelerator experiments 
probing the — > v e appearance channel, which is governed by the (Am 2 , #13) parameters (although with an additional 
dependence on #23 and S, absent in short-baseline reactor experiments). In particular, the Tokai-to-Kamioka (T2K) 
experiment has observed 6 electron-like events with an estimated background of 1.5 events, rejecting #13 = at the 
level of 2.5(7 [23[. The low background level makes the T2K results particularly important and robust. Shortly after, 
the Main Injector Neutrino Oscillation Search (MINOS) experiment has reported the observation of 62 electron-like 
events with an estimated background of 49 events, disfavoring 6*13 = at 1.5cr [2^, [2f|. Taken together, these data 
suggest sin 2 #13 ~ few %, in agreement with our previous hints Q discussed above. 

It makes then sense to update our previous global analyses of oscillation data 0, H, [ll| by including the latest T2K 
and MINOS results, as well as other data which have been published in the last few years. Our main result is 

, f 0.021 ± 0.007 , old reactor fluxes 

sin 2 0i3 = <^ (la) , 3) 

\ 0.025 ±0.007, new reactor fluxes y J J 

corresponding to a > 3er evidence in favor of nonzero $13 (while previous hints did not exceed the ~ 2er level). We 
discuss below some details of our current analysis, and a few relevant implications of $13 > for neutrino physics. 
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II. INPUT DATA AND THEIR ANALYSIS 

We briefly summarize the main updates in our global analysis of neutrino oscillations driven by Sm 2 and by Am 2 , 
respectively, with respect to our previous works on the subject. 

A. Data sensitive to 5m 2 

These data include events collected in solar neutrino experiments and in the long-baseline reactor experiment 
KamLAND. As described in [|| [26|, we take as free parameters (5m 2 , 612, 6*13), as well as the four geoneutrino event 
rates corresponding to decays of natural Thorium and Uranium in the KamLAND and Borexino experiments, which 
are then marginalized away. With respect to (2(| , wc update the KamLAND results on reactor and geoneutrino data 
[3113] , and include the Borexino total event rates of 7 Be [28| and 8 B [29[ solar neutrinos. As a side result, we obtain 
good agreement with the latest geoneutrino analysis performed by the KamLAND collaboration [27| ■ 

The interpretation of both solar and reactor neutrino data depend, to some extent, on the estimated fluxes (and 
their uncertainties) at their sources, namely, the solar and reactor cores, respectively. Concerning solar neutrinos, 
the long-standing issue of high versus low metallicity (see [3(| for a review of the subject) affects only marginally the 
estimates of the oscillation parameters Hp . For instance, we find that by using the standard solar model (SSM) 
inputs GS98 and AGSS09 from Table 3 of [30|, the significance of the solar+KamLAND constraints on sin 2 #13 changes 
by only ~ O.lcr. For definiteness, we use the reference BP05(OP) SSM from Table 2 of [31], which gives fit results 
intermediate between the GS98 or AGSS09 SSM. 

Conversely, a more relevant issue has been recently pointed out in the context of reactor neutrino fluxes, whose 
detailed reevaluation in [l8| suggests an average increase of about 3.5% in normalization, with respect to previously 
accepted standards [32| (see also [I(|). The increase would then indicate extra electron flavor disappearance in 
KamLAND [l9| which, in the context of Zv oscillations, are expected to induce a small but nonnegligible shift in #12 
and 0i3 [22j. We will thus present results obtained by using either "old" [32| and "new" [l8[ reactor fluxes and their 
uncertainties in the analysis of reactor data, in order to illustrate the effect of such flux systematics. 

B. Data sensitive to Am 2 

These data include, in our analysis, events collected in the atmospheric neutrino experiment Super-Kamiokandc 
(SK), in the short-baseline reactor experiment CHOOZ, and in the long-baseline experiments K2K, T2K, and MINOS. 

With respect to our previous analysis Q, the CHOOZ data [IH are analyzed also with new input reactor fluxes 
fl8j , whose higher normalization leads to a small disappearance effects which slightly favors nonzero #13 and weakens 
its upper bounds fl9j |. 

We update the atmospheric neutrino analysis @, Q by including SK-II and SK-III data [l4], [13] , together with 
improved estimates for the associated systematic errors [3J]. However, a reduction of the SK-I+II+III data is un- 
avoidable, since the official SK analysis includes many categories of events, bins, and systematics, which cannot be fully 
reproduced outside the collaboration, as already noticed [22[. In particular, we exclude from our analysis multi-ring 
and pion-decay events, and we group together some classes of partially contained events (stopping and through-going) 
and of upgoing muon events (showering and nonshowering) . We think that, despite some unavoidable approximations 
and loss of information, our analysis remains a useful, independent study of SK atmospheric data within a full Zv 



oscillation framework. A similar attitude has been adopted in [13j, while the authors in [21[ have directly taken the 
official SK likelihood function for the (Am 2 , #23, #13) parameters from a 2>v analysis with 5m 2 set to zero 14[. As a 
consequence, the atmospheric neutrino likelihood adopted in (2l| lacks the 3v sublcading oscillation effects driven by 
(5m 2 [iH , which are potentially relevant at the current level of accuracy 0, G3] ■ 

Concerning long-baseline accelerator data in the — > disappearance channel, we maintain our previous K2K 
analysis without changes @, @], but we update the MINOS spectrum analysis from (36|. We do not include MINOS 
antincutrino disappearance data, which seem to indicate a puzzling deviation from the neutrino best-fit parameters 
(Am 2 , #23 ) [37} • The deviation should be carefully monitored in the future, since it cannot be interpreted in the 
standard iv framework. 

In the appearance channel — > v e , we include the recent results from [2~il [25| . by fitting the total (spectrum- 
integrated) rate of electron-like events at the far detector, taking into account its statistical and systematic errors. 
We reproduce with good accuracy the marginalized MINOS confidence level contours (25| in the plane charted by the 
CP phase 8 and by the dominant mixing parameter sin 2 #23 sin 2 2#i3, for both normal and inverted hierarchy. Last, 
but definitely not least, the crucial T2K appearance data recently reported in [23| are included by fitting the total 
(spectrum integrated) rate at the far detector within Poisson statistics. We reproduce with good accuracy the T2K 
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confidence level contours in the plane charted by the CP phase 5 and by the mixing parameter sin 2 2#i3 for fixed 
values of the other oscillation parameters, in both normal and inverted hierarchy (23[. [Note that all such fixings are 
removed in the global analysis below.] 

We emphasize that the official appearance data analyses in T2K [23| and MINOS [25[ include also spectral infor- 
mation and uncertainties, which are quite difficult to model outside the respective collaborations. However, we think 
that their current parameter fits can be reproduced, to a large extent, by a single number: the total appearance rate 
(normalized to published no-oscillation expectations in each experiment). In particular, we have not found appre- 
ciable variations by performing trial fits with binned spectra. We also note that the continuous degeneracy of the 
(&13,5) minima, trivially expected in a total-rate fit, is not lifted by the additional spectral information of complete 
fits [23, Hil- Of course, future, higher-statistics T2K and MINOS results might entail more constraining spectral 
information than in current data. 

C. Three- neutrino analysis 

In (5m 2 -sensitive oscillation searches, the relevant 3v variables are (5m 2 , sin 2 612, sin 2 #13). A very minor dependence 
of solar neutrinos on ±Ato 2 0, [38| is also kept for the sake of precision. 

In Am 2 -sensitive oscillation searches, we take as free parameters (±Am 2 , sin 2 #23, sin 2 #13), while (5m 2 , #12) are 
fixed at their best-fit values from the analysis of <5m 2 -sensitive data. Our atmospheric v analysis is restricted to the 
two CP-conserving cases cos 5 = ±1 @, @|; we plan to remove such restriction in the future, so as to analyze both 
atmospheric and long-baseline accelerator data with unconstrained 5 (see [l3[ for a previous analysis of this kind). 
Anyway, as argued at the end of Sec. Ill, the evidence of #13 > can only be strengthened for unconstrained 5. 

In the global fit, note that the common parameter sin 2 #13 is constrained by both classes of oscillation searches. As 
in previous works 0, [|| fill ] , we will present bounds on the continuous parameters 

(5m 2 , sin 2 8 12 , sin 2 13 , sin 2 6> 23 , Am 2 ) , (4) 

marginalized over the four discrete cases 

[sign(Am 2 ) = ±1] <g) [cos 5 = ±1] , (5) 

unless otherwise stated. Concerning Am 2 , note that our convention in Eq. ^ absorbs the trivial 5m 2 difference 
arising in the best-fit values (for normal and inverted hierarchy) of the alternative conventional mass parameters 
Am^ = m\ — m\ or Am 2 2 = m 2 — m\. 

A final remark is in order, concerning some issues recently raised by a reevaluation of reactor fluxes [HI, [l9j]. The 
higher normalization suggested in [l8j produces extra disappearance in KamLAND and CHOOZ, which jointly favor 
higher #13 and #12 in a 3v framework. However, it also produces a small disappearance effect in very short-baseline 
experiments (not included in our fit) which is at variance with standard 3v oscillations, and suggests hypothetical 
sterile neutrino oscillation driven by a mass gap AM 2 ~ 0(1) cV 2 and by relatively small mixing with v e [l9|, |39| — 
H2|. Conversely, the old normalization is essentially consistent with the 3v scenario, but slightly weakens the global 
indications for 6*13 > 0. Our educated guess (or prejudice) is that the true normalization may be intermediate between 
the old and the new one, and that previous normalization errors were underestimated. Anyway, we shall show results 
for both old and new reactor fluxes, their difference being an indication of the corresponding systematic uncertainties. 
We emphasize that the understanding of these puzzling systematics definitely requires further experimental data. We 
also note that fluxes even lower than the "old" ones have been investigated as a possible phenomenological explanation 
of very short baseline reactor data 0, In general, the lower the reactor fluxes, the lower the allowed values of 
#13, as discussed below. 

III. RESULTS 

Figure 1 shows the main results of our global fit, in terms of allowed ranges for each of the oscillation parameters 
in Eq. (j4|), as obtained by marginalizing the remaining parameters over all the four discrete cases in Eq. (0. The 
vertical scale represents the number of standard deviations (Na = y 7 A% 2 ) from the best fit point. The lines departing 
from each best fit would be perfectly straight and symmetric for exactly gaussian error distributions. This is nearly 
the case for all but the sin 2 623 parameter, which shows a skewed preference for the first octant at the level of la, in 
qualitative agreement with our previous analyses 0, The estimates of sin 2 #13 and sin 2 0i 2 are affected by reactor 
flux systematics; in particular, the dashed lines refer to the analysis with new reactor fluxes [l8j . which turn out to 
shift both parameters by roughly +0.005 with respect to the case with old reactor fluxes (solid lines). In both cases, 
however, there is a clear evidence in favor of sin 2 6*13 > 0, at a confidence level of at least 3a. This evidence is the 
most important outcome of our work. 
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Synopsis of global 3v oscillation analysis 




7.0 7.5 8.0 0.25 0.30 0.35 0.00 0.02 0.04 0.3 0.4 0.5 0.6 0.7 2.0 2.5 3.0 



5m 2 /10~ 5 eV 2 sin 2 9 12 sir|2e i3 sin2 e 23 Am 2 /10" 3 eV 2 

FIG. 1: Global 3v oscillation analysis. Bounds on the mass-mixing oscillation parameters, in terms of standard deviations from 
the best fit. Solid (dashed) lines refer to the case with old (new) reactor neutrino fluxes. Note the > 3cr evidence for #13 > 0. 

Figure 2 breaks down the global evidence for sin 2 #13 > into two separate contributions coming from the data 
sets sensitive to either 5m 2 (Solar+KamLAND) or Am 2 (ATM+LBL+CHOOZ), assuming old and new reactor fluxes 
(left and right panels, respectively). Remarkably, the two data sets agree very well, with best fits rather close to each 
other in both panels, and with nearly gaussian uncertainties in all cases. The bounds from combined (ALL) data 
appear to be currently dominated by Am 2 -sensitive experiments — not surprisingly, since the T2K appearance results 
alone account for more than 2a }23| . The T2K experiment, currently limited by statistics rather than by systematics, 
is expected to improve significantly the bounds on 813 in future physics runs [23| . We also find it useful to summarize 
the ±ler ranges of sin 2 6*13 in a different format in Fig. 3, where the solid and dashed error bars refer to old and new 
reactor neutrino fluxes, respectively. 



Global evidence for 6 13 > 




FIG. 2: Breakdown of the evidence for #13 > from the global fit (ALL) into contributions coming from 5m 2 -sensitive data 
(Solar+KamLAND) and from Am 2 -sensitive data (ATM+LBL+CHOOZ). The left and right panels refer to old and new fluxes, 
respectively. 
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Global evidence for 9 13 > 
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FIG. 3: Global 3v analysis. Preferred ±1<t ranges for the mixing parameter sin 2 #13 from partial and global data sets. Solid 
and dashed error bars refer to old and new reactor neutrino fluxes, respectively. 

Tabic I reports the bounds shown in Figs. 1-3 in numerical form. All the bounds arc largely uncorrclatcd from each 
other; e.g., the allowed ranges of Sm 2 and Am 2 are basically independent on variations of the mixing angles within 
their uncertainties (not shown). Nevertheless, we find it useful to report the joint ranges for the mixing parameters 
sin 2 8ij . which can be used to test specific predictions of theoretical models for neutrino mixing, and which allow to 
highlight the impact of recent appearance data. 

Figure 4 shows the joint contours at 1, 2 and 3ct (A% 2 = 1,4 and 9) for each possible couple of sin 2 Oij parameters, 
in the analysis with old reactor fluxes. Including new fluxes, the best fits and the associated Na contours are all 
translated by small amounts (< la) indicated by arrows. As a result of the dominance of T2K data in the 9 13 fit, the 
correlation in the (sin 2 # 12 , sin 2 13 ) plane induced by (5m 2 -sensitive data [3, [TTJ] is no longer apparent. Conversely, 
there is a weak anticorrelation in the (sin 2 #23, sin 2 6*13) plane for relatively high #13, due to the fact that the long- 
baseline — > v e appearance probability is dominated by the product \U^U e3 \ 2 oc sin 2 #23 sin 2 2#i 3 . If we had we 
hypothetically obtained a best fit value sin 2 #23 — 0.5 instead of 0.42, the best-fit value of sin 2 6*13 would have been 
presumably slightly lower (by about —0.002, according to our educated guess). 



TABLE I: Results of the global 3v oscillation analysis, in terms of best-fit values and allowed 1, 2 and 3a ranges for the 
mass-mixing parameters, assuming old reactor neutrino fluxes. By using new reactor fluxes, the corresponding best fits and 
ranges for sin 2 812 and sin 2 #13 (in parentheses) are basically shifted by about +0.006 and +0.004, respectively, while the other 
parameters are essentially unchanged. All the results are marginalized over the four discrete cases in Eq. (0). 



Parameter 


<5m 2 /10- 5 eV 2 


sin 2 612 


sin 2 #13 


sin 2 623 


Am 2 /10 -3 eV 2 


Best fit 


7.58 


0.306 


0.021 


0.42 


2.35 






(0.312) 


(0.025) 






lcr range 


7.32 - 7.80 


0.291 - 0.324 


0.013 - 0.028 


0.39 - 0.50 


2.26 - 2.47 






(0.296 - 0.329) 


(0.018 - 0.032) 






2cr range 


7.16 - 7.99 


0.275 - 0.342 


0.008 - 0.036 


0.36 - 0.60 


2.17 - 2.57 






(0.280 - 0.347) 


(0.012 - 0.041) 






3(7 range 


6.99 - 8.18 


0.259 - 0.359 


0.001 - 0.044 


0.34 - 0.64 


2.06 - 2.67 






(0.265 - 0.364) 


(0.005 - 0.050) 
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sin 2 12 sin 2 23 

FIG. 4: Global 3v analysis: Joint contours at 1, 2 and 3a (A% 2 = 1,4 and 9) for couples of sin 2 9ij parameters, assuming old 
reactor neutrino fluxes. For new reactor fluxes, the best fits (and, to a large extent, also the contours) are shifted as indicated 
by the arrows. 



All the results reported in Figs. 1-4 and in Table I are marginalized over the four discrete cases in Eq. ([5]). For 
completeness, we also show in Fig. 5 the breakdown of the #13 bounds over these four options (for the case of old 
reactor fluxes only), in terms of standard deviation ranges for the parameter cos 5 sin #13 = ± sin #13 in both normal 
and inverted hierarchy, for the relevant data set ATM+LBL+CHOOZ (while <5m 2 -sensitive data are ^-independent). 
We do not find any significant difference in the position and likelihood of the best fits points for normal and inverted 
hierarchy. However, we find, as in our previous analysis 0], a preference for the CP-conserving case cos 5 = — 1 versus 
+1. In [2|] we argued that the weak preference for both sin 2 6*13 > and cosi5 < were tied by the interference term 
driven by 5m 2 and #13 [35lj in the atmospheric neutrino oscillation probability. The preference for cos 6 sin #13 < 
in Fig. 5 appears to be more pronounced (~ 2er) than in 0, presumably because it correlates with the (now more 
robust) preference for sin 2 #13 > 0. It is intriguing to notice that a weak preference for cos S sin 6*13 < has also been 
found in the preliminary atmospheric 3i> analysis from the SK collaboration (l7j . Time will tell if this is another hint 
or mere fluctuation. 



ATM + LBL + CHOOZ (Am 2 , 8 23 free) 
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FIG. 5: Standard deviations from the best fit in terms of the variable cos <5 sin #13 for the two CP parities (cos<5 = ±1) and 
for both normal hierarchy (NH) and inverted hierarchy (IH), using the ATM+LBL+CHOOZ data set. The curves refer to the 
analysis with old reactor neutrino fluxes; similar results (not shown) are obtained for new reactor fluxes. 
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As a final remark, we note that the evidence in favor of 6*13 > 0, found at a level > 3er in this paper, could only 
be strengthened by taking the CP phase S as a free, continuous parameter. Indeed, for #13 = 0, the value of 5 is 
physically irrelevant, and thus the value of Xo = X 2 |ei 3 =o remains unchanged. For #13 > and unconstrained 5, the 
absolute minimum of x 2 can only decrease (as compared to our constrained cases, cos S = ±1), leading to an increase 
of Ax 2 = Xo ~ Xmim an d thus to an even higher statistical evidence in favor of 6*13 > 0. Of course, for free 6 the best 
fit of sin 2 #13 might be reached at values slightly different than those discussed herein (0.021 or 0.025, depending on 
reactor systematics). 

IV. DISCUSSION AND CONCLUSIONS 

The hints of 613 > that wc pointed out in [3, 0] are consistent with recent T2K (H] and MINOS [H] data 
and, in combination, provide an evidence for $13 > at the level of > 3o\ Such evidence, and the preference for 
values in the range sin 2 #13 ~ 0.01-0.04 at 2<r (see Table I) may have far-reaching consequences in neutrino physics. 
First of all, it is crucial to experimentally test these findings, not only with further long-baseline appearance data at 
accelerators, but also with with short- baseline disappearance searches at reactors [l5j |. If confirmed, the evidence for 
sin 2 9 13 ~ few % would open the door to CP violation searches in the neutrino sector, with profound implications for our 
understanding of the matter-antimatter asymmetry in the universe. At the same time, matter effects on (±Am 2 , #13)- 
driven oscillations appear now as a much more promising tool, from both the experimental and theoretical viewpoint, 
to derive indications about the hierarchy from neutrino propagation in the Earth or in supernovae. Concerning 
nonoscillation searches, the evidence for #13 > provides small but "guaranteed" contributions of the third neutrino 
mass ?7i3 to both single beta and double beta decay searches, which need to be accounted for in detailed analyses. 
From a more theoretical viewpoint, relatively large values for #13 will certainly trigger new ideas for model building. 
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